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SIMULATION OF HYPERSONIC SCRAMJET EXHAUST 


By Ro A. Oman, K. M„ Foreman, J. Leng, and H. B. Hopkins 

Research Department 
Grumman Aerospace Corporation 
Bethpage, New York 11714 


INTRODUCTION AND SUMMARY 


This report describes a plan and some preliminary analysis 
for the accurate simulation of pressure distributions on the 
afterbody/nozzle portions of a hypersonic scramjet vehicle such 
as shown in the figure on page 2. It is apparent from the cur- 
rent design philosophy for these vehicles that the scramjet ex- 
haust gases will play a major role in the determination of the 
aerodynamic forces and moments acting on the vehicle. Conven- 
tional aerodynamic testing methods offer no straightforward ap- 
proach that will meet all of the vehicle development require- 
ments. In this first investigation our objectives were to es- 
tablish the standards of similitude for a hydrogen/air scramjet 
exhaust interacting with a vehicle afterbody, determine an ex- 
perimental technique for validation of the procedures that will 
be used in conventional wind tunnel facilities, suggest a pro- 
gram of experiments for proof of the concept, and explore any un- 
resolved problems in the proposed simulation procedures. We have 
fulfilled all of these objectives, and the second phase activi- 
ties of proof-of -concept experiments and more detailed analysis 
are ready to begin. 

The over-all plan for vehicle development that we have de- 
fined during the program covered by this report has four sepa- 
rate activities that culminate in a final determination of the 
afterbody pressure distribution in a wind tunnel. First, the 
expansion characteristics of the hydrogen/air combustion products 
at true flight enthalpy will be determined by pressure and other 
diagnostic measurements in a special facility called a detonation 
tube simulator. Second, a small number of substitute gases will 
be tested in the same facility to show that they are capable of 
reproducing at moderate temperature levels the pressure distribu- 
tions of the hot combustion products. Third, a series of experi- 
ments will be conducted with substitute gases and with combustion 
products to develop methods for assessing such nonuniformities as 
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combustor partition wakes, combustor wall boundary layers, cor- 
ner flows, and combustor end effects. Fourth, the qualified 
substitute gases will be used in vehicle and combustor develop- 
ment testing in conventional wind tunnel facilities. The four 
activities can, in part, run concurrently. 

One of the most difficult simulation problems for the 
scramjet exhaust is to reproduce the effect of the correct chem- 
istry along with the high total enthalpy (flight compression 
plus combustion heating) in the exhaust flow. To achieve this 
combination with the true combustion mixture in a conventional 
wind tunnel model of reasonable size and complexity is a virtual 
impossibility. The approach we have taken in the current work 
is to investigate the use of a cold substitute gas, or blend of 
gases, to simulate the scramjet exhaust. In order to select 
such a substitute gas, we must have a comparison technique for 
validating the gas mixtures chosen to represent the exhaust at 
different flight conditions. The best validation will come from 
measurement of actual pressure distributions over a scale model 
of the vehicle afterbody with an exhaust gas consisting of the 
products of combustion of hydrogen and air at an enthalpy level 
as close to that of the real engine as possible. These pressure 
distributions can then be compared to those obtained with selec- 
ted substitute gas mixtures. 

An extensive analysis of the requirements that must be ful- 
filled to ensure fidelity of model pressure distributions to the 
flight case showed that the most critical parameters were match- 
ing combustor exit Mach number and following the variation of 
specific heat ratio with increasing M in an isentropic expan- 
sion. Many other important simulation requirements have also 
been satisfied by the plan we have devised. The requirements of 
pressure simulation within the diverging contact surfaces of the 
underexpanded exhaust are met to a remarkable degree. Reynolds 
number is matched, Prandtl number is quite closely matched, 
and the viscosity variation through the boundary layer is shown 
to be nearly that needed for an ideal boundary layer match. 
Slight adjustments in the substitute gas blends can easily be 
made to allow for different degrees of chemical equilibrium in 
the flight exhaust. Even the strict similitude requirements of 
the turbulent mixing layer between exhaust and external flow 
could be fulfilled by the substitute gases if they were used in 
a special Freon wind tunnel — but we recommend a separate in- 
vestigation of hypersonic mixing layer similitude for the pur- 
pose of easing this requirement. In short, the substitute gas 



concept looks like a valid approach, and the next step is to 
prove it by comparison to a valid standard flow. 

The only known technique for producipg the correct enthalpy 
and chemistry in the laboratory appears to be the detonation 
tube concept pioneered by Grumman for the LM program, and used 
subsequently to simulate the exhaust from the hydrogen/oxygen 
main engines for the Space Shuttle (Ref. 1). This technique is 
ideally suited to scramjet simulation because it combines shock 
heating of the gas with the heat of combustion to obtain the re- 
quired enthalpy. One of the main objectives of our current work 
was to investigate analytically the operating range of the deto- 
nation tube facility and to match it to the proposed flight en- 
velope of the scramjet aircraft. We have done this and shown 
the technique to be capable of high quality simulation through- 
out the complete flight regime and beyond. Having shown that 
operating range is not a significant problem for the detonation 
tube system, we selected a nominal Mach 8 flight condition as a 
baseline condition for the first set of demonstration tests. 
Picking a single flight Mach number restricts the required model 
building to a set with one exit area ratio yet allows a full 
range of experiments on pressure and enthalpy levels, substitute 
gas matching, and exit plane partition wakes. 

The next objective was to examine and select possible sub- 
stitute gases for subsequent experiments. Three promising bi- 
nary mixtures were chosen, and two dimensional pressure distri- 
butions were calculated for the flight baseline conditions on 
the real engine, the detonation tube simulation, and the three 
substitute gases. The selection procedures involved a one di- 
mensional screening technique, after which the two dimensional 
method of characteristics calculations were made to compare ac- 
tual pressure distributions over the two dimensional representa- 
tion of the geometry that we propose to test in the next phase. 
On page 5, we show a preview of Fig. 16, which describes better 
than words the matches of axial and normal forces that we can 
get. The blends used are easily prepared and modified, safe, 
and relatively inexpensive. Although only inviscid flows are 
shown, the friction drag will also be properly scaled when nor- 
malized by the exit plane pressure. 

In summary, the results of the present investigation show 
that the detonation tube simulator should provide an excellent 
simulation of the hydrogen/air scramjet exhaust flow over the 
vehicle afterbody for the entire flight regime. There are 



several nontoxic , relatively inexpensive substitute gas mixtures 
available for use in wind tunnels that will give adequate repre- 
senations of the pressure distributions over the afterbody due 
to the scramjet exhaust. In order to choose the best substitute 
gas mixture for a given flight regime, comparisons with the de- 
tonation tube results and with mixing layer similitude require- 
ments will be made. An experimental program to prove this con- 
cept at a representative point in the flight envelope is pre- 
sented. We expect that such a program will provide the ability 
to simulate correctly the aerodynamic loads due to the engine 
exhaust in subsequent experiments in Langley facilities. 
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SYMBOLS 


A 


Ar 


D 


12 


D, 


D III 

Ec 


Freon 12 
Freon 13 
Freon 13B1 
Freon 14 
Freon E5 


H 

HFA 

K 

k 

L 

LM 


area 

argon atom 
sound speed 

specific heat at constant pressure 

specific heat at constant volume 

binary diffusion coefficient, defined on page 64 

first Damkohler parameter, defined on page 15 

third Damkohler parameter, defined on page 16 

Eckert number, defined on page 16 

fluorine atom; also abreviation for Freon 

dichlorodifluoromethane (CCI 2 F 2 ) 

chlorotrif luorome thane (CCIF^) 

bromotrifluoromethane (CBrF^) 

tetraf luorome thane (CF^) 

a complex fluorinated ether 

enthalpy; also hydrogen atom 

hydrogen fueled aircraft 

Kelvin temperature scale 

thermal conductivity 

length 

NASA/Grumman Lunar Module 



Le 

M 

% 

m 

N 

Nu 

n 

0 

P 

Pr 

• 

q 

R 

Re 

S 

SCRJ 

Sc 

T 

t 

U 

V 

X 

X 


Lewis number, defined on page 15 
Mach number (= V/a) 
molecular weight 
unit of length (meter) 

nitrogen atom; also unit of force (Newton) 
Nusselt number, defined on page 17 
exponent in temperature-viscosity power law 
oxygen atom 
pressure 

Prandtl number, defined on page 1A 
heat transfer rate 
universal gas constant 
Reynolds number (= pVL/|i) 
sulphur atom; also entropy 
Supersonic Combustion Ramjet 
Schmidt number, defined on page 15 
temperature 
time 

velocity (in detonation tube) 
velocity 

coefficient in chemical mixture equation 

(pp. 34 and 61); also afterbody axial distance 

nondimens ional distance (= X/Y^) 
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y 

Y 

a 

P 

y 

e 

[!■ 

p 
a . 

T 

<P 


coefficient in chemical mixture equation, page 34; 
also afterbody radial distance 

nondimens ional distance (= Y/Y^) 

angle of attack 

inlet ramp angle 

ratio of specific heats, c / c 

average y between two frozen thermodynamic states 

average y between two equilibrium thermodynamic 
states 

temperature difference between two points 
viscosity, or prefix "micro" (10 
density 

thermodynamic function defined on page 34 ; also 
molecular collision diameter constant defined on 
page 64 

ratio of two different absolute temperatures 
equivalence ratio 

mixture viscosity parameter defined on page 62 


Subscripts 


o 


refers to total (stagnation) condition 


1 


refers to an arbitrary thermodynamic state point 


2 / refers to an arbitrary thermodynamic state point; 

also thermodynamic state behind incident shock wave 
in detonation tube 
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refers to combustor exit plane 

refers to unexpanded, high pressure driver gas 

refers . to stagnation condition or region behind re- 
flected shock wave in detonation tube 

refers to undisturbed free stream 

refers to Chapman-Jouguet thermodynamic state 
th 

refers to i chemical specie; also internal flow 

refers to external flow; also shock tube nozzle 
exit plane 

refers to a mixture of gases 
refers to constant entropy 


Superscript 


refers to sonic condition 
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SIMILITUDE AND SUBSTITUTE GASES 


The purpose of this part of the investigation is to identi- 
fy the thermodynamic and transport property relationships neces- 
sary to ensure that the flow fields of scaled test nozzles will 
produce the same pressure distributions as those of the full 
scale exhaust nozzles used in supersonic combustion ramjets 
(SCRJ). Because of the method selected to achieve this result, 
a critical associated task is the selection of suitable test 
gases to be substituted in wind tunnel tests for the extremely 
high temperature exhaust gases produced by the hydrogen fueled 
aircraft (HFA) . 

The nondimensional similarity parameters for inviscid and 
viscous flow are computed for the SCRJ and used to formulate 
substitute gases with the required properties, assuming geomet- 
ric similarity between the model and full scale nozzle configu- 
rations. A significant effort has been required to compile the 
most current and authoritative transport properties for the 
likely substitute gas constituents and to determine the most 
accurate combining rules for estimating transport properties of 
the gas mixtures. Sixteen gases considered safe, easily han- 
dled, chemically stable, relatively inexpensive, and commercial- 
ly available were considered. The most promising substitute gas 
blends (binary mixtures) from those sixteen were determined by 
the following four-step matching process. Initially, for spe- 
cific flight conditions, the combustor exit specific heat ratio, 
73 , and the approximate range of 7 variation in the SCRJ ex- 
haust gas expansion are used to screen the candidates. Then, 
the transport similarity parameters (such as Reynolds, Prandtl, 
and Schmidt numbers) are used to further narrow the choices of 
substitute gas components. Subsequently, the specific initial 
temperature and pressure conditions and fractional composition 
of the substitute gases that exhibit the best pressure distri- 
bution match to the SCRJ gases along the entire nozzle expansion 
process are selected, assuming one dimensional isentropic flow 
and the particular needs or restrictions of the testing tech- 
nique to be used. Finally, the best matches of the one dimen- 
sional computations are confirmed by two dimensional analysis. 

As a result of this procedure we have concluded that binary 
mixtures of four Freon compounds (Freon 12, 13, 13B1, and 14) 
with argon additive are good substitutes for SCRJ combustion 
gases for facility testing of model exhaust nozzles at condi- 
tions simulating flight between Mach 4 and 10. Using these 
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mixtures with about £ scale models provides complete dynamic 
inviscid simulation of the flow field pressure distribution 
along the SCRJ nozzle walls and will closely approximate trans- 
port similitude as well. These mixtures are compatible with 
either the wind tunnel or the detonation tube test facilities. 

Similitude Parameters 

The nondimens ional parameter requirements for dynamic simi- 
larity of models in fluid dynamic systems without heat transfer 
and chemical reactions have been well developed in the litera- 
ture (e.g., Ref. 2). For inviscid, two or three dimensional com- 
pressible flow and geometrically similar systems, the Mach number 
of the model and prototype must .be reproduced throughout the 
flows so that geometric relationships between physical surfaces 
and the characteristic directions of propagation in the flow are 
preserved. A primary additional requirement for this geometric 
preservation in a nozzle exit region is that the ratio of pres- 
sures between the external and internal flows is established such 
that all wave systems that affect regions of measurement are 
properly reproduced. In most cases this means a simple matching 
of pressure ratio, but in some situations corrections for gas 
thermodynamics are required. If the region of interest is re- 
stricted to that inside the Mach cones from the exit lips and 
therej are no compression waves from the exit lips, it is neces- 
sary to ensure only that the external pressure is lower than the 
corresponding flight external pressure. For perfect gases in 
equilibrium it is then also necessary and sufficient for inviscid 
simulation that the ratio of specific heats, y, be the same in 
model test and in prototype flight. However, the gases involved 
in SCRJ flows show very significant variations in specific heat 
within the flow, but in the exhaust region of interest they are 
well represented by the equation of state P= pRT/^, where R 
is the universal gas constant and the molecular weight #1 can be 
regarded as a local constant (see additional discussion below). 
This type of behavior is referred to as thermally perfect, and 
our subsequent discussion and analysis will employ that simpli- 
fication. , 

We will show below that a thermally perfect gas flow will 
be dynamically similar to another flow for matched geometries as 
defined above whenever the local ratio of specific heats, y , 
is the same at corresponding geometric points in the two flows 
and molecular weight changes along the streamlines of each flow 
are negligible. Because temperature levels in the two flows 



(model and prototype) may be radically different, the proper way 
to ensure this match is to ensure that 7(T/T re f) is the same 
function for both gases, where Tref is the temperature at any 
convenient geometric reference point common to the two flows. 

We have chosen the combustor exit plane temperature, T3 , as 
our reference temperature. The y function used here is the 
local ratio of specific heats, c p /c v , which should not be 
confused with the isentropic exponent or the effective y that 
can be used to characterize a shock or other finite property 
change of a process. These latter effective 7' s characterize 
a process in terms of an equivalent process for which the local 
specific heat ratio would be unchanging, while the 7 we use 
changes throughout the flows. We will henceforth use y only 
to characterize relationships among gas properties at a point 
or for infinitesimal processes, not across finite state changes. 
For finite processes we will use 7. 

The flows of greatest concern in SCRJ simulation are either 
locally isentropic, or include shock processes. Starting with 
the shock case, we can write the Rankine-Hugoniot equations to 
relate properties on either side of a shock wave 

p V = p,V, 

•a a K b b 


P a + »£ - P b + pX 
"a + i ^ - «b + i i 


Normalizing all properties by their conditions at the exit 
plane (subscript 3 ), the local velocity can be expressed in 
terms of the local Mach number 


J 1 x 


V = V-R X 3 V - M 


and the total enthalpy is expressed 


H 


c dT = RT,, 
P 3 


Y 


ftOv-i) 


<8 
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The above system of relationships is closed by the equa- 
tion of state 



so that any property ratio at b is determinable by knowledge 
of the state a and the functions y(T/T 3 ) and #!(T/T 3 ). Note 
that although fn. is not a function of temperature ratio alone 
it is in a particular flow system where p and T are coupled 
to a particular relationship. 


Similarly, for external flow the continuity constraint 
equation is set by the geometry and the constraint of constant 
entropy is added 



0 


The internal energy E can be written 


E = 


c v dT = RT 3 


1 J_ TN 

57!Oy- 1) a \T,/‘ 


so we get 


P 3 ra 3 


*1(7-1) 


it) 




= constant 


Thus we see that if the Mach number at the boundaries is 
everywhere the same, matching of the functions y(T/T 3 ) and 
#;(T/T 3 ) is sufficient to ensure that all property ratios 
(e.g., P/P 3 , p/p 3 , etc.) are matched in two different inviscid 
flows, with or without shock processes. 


For the M^ = 10 flight case we find that with no shock 
waves present, the equilibrium molecular weight varies from 
about 24.3 at the combustor exit to 24.6 at the end of the 
afterbody. Even if the flow were shocked back to the sonic con- 
dition, the molecular weight would not go below 22.7. For the 
Moo = 8 flight case, the variations are even less. We have 



therefore concluded that the molecular weight can be treated as 
a constant in the SCRJ exhaust, even in the presence of shocks, 
and that substitute gases need only have matching y and Mach 
number distributions for inviscid similarity. Note again that 
77i of the substitute gas and SCRJ flow need not be the same. 

Inclusion of viscous regions in the similitude requires 
matching of Reynolds number, Re. Also necessary for dynamic 
similarity in either the viscid or inviscid case is the similari- 
ty of properties of the gas used for model testing and the gas 
associated with the prototype. Principal inviscid gas proper- 
ties are represented by Prandtl number, Pr, and the variation 
of viscosity and thermal conductivity with temperature, where 
Pr = Cp M-k”!, which can be considered the ratio of momentum 
diffusion to thermal diffusion. 

The viscosity of gases is independent of pressure at rela- 
tively low pressure and generally increases with temperature. 

For nonpolar gas molecules, the viscosity varies as a power of 
temperature, n, where n approaches 0.5 for high tempera- 
ture gases and is nearer 1.0 for low temperatures. Polar 
gases, such as the H 2 O which comprises approximately one third 
of the SCRJ exhaust gas composition, cannot have their viscosity 
variation so easily characterized. However, over a finite tem- 
perature range of interest, the viscosity calculations or data 
for SCRJ and test gases can be fitted locally by a law such as 
H-i/m- 2 = (Ti/T 2 ) n > or by a two-constant Sutherland equation such 
as |i = bT^/^^s + T)~ , where b and S are constants. In 
this way the model test and prototype operation conditions can 
be examined for similarity of the viscosity distributions. 

Thermal conductivity of gases, k, increases nearly linearly 
with temperature. However, for diatomic and polyatomic gases 
no general temperature relation has been found satisfactory. 

Curve fitting of k to a polynomial in temperature using test 
or computed data is one way of representing this transport pro- 
perty. The use of Pr as a similarity parameter characterizes 
the relationship between 1 1 and k, but the relationship is 
strictly correct only if Pr is invariant with temperature. 

In combustion, chemical reactions, and heat transfer pro- 
cesses, four additional independent similarity parameters are 
necessary. The diffusion of mass from one part of the multi- 
specie system to another must be considered. The Schmidt num- 
ber, Sc, ratios the diffusion of momentum to mass 



From Sc and Pr the Lewis number, Le, is derived to give the 
ratio of mass diffusion to thermal diffusion: 

Le = = P^c D, 0 (RTk)" 1 . 

Sc p 12' 

The first and third Damkohler parameters, Dj- and Djjj 
also are important, especially in chemically reacting flows. The 
ratio of residence time of a reacting gas mixture in a flow chan- 
nel, L/V, to the duration of the controlling chemical reaction, 
t^, is given by the first Damkohler grouping, Bj = L/Vt^. For 
Dl » 1, the gas can be considered in local equilibrium even 
though several secondary reaction products may not have reached 
final concentrations. For Dj « 1 the reaction products will 
be nearly frozen in composition during the process. For accurate 
details of the flow chemical kinetic rate data must be used with- 
the flow equation. In both the SCRJ exhaust gas flow and detona- 
tion tube simulator, the kinetics can be important in determining 
the flow field and pressure distribution for high Mach number 
(altitude) cases, although the coupling to pressure is rather 
weak. For the binary substitute gas mixture, the components are 
chemically unchanging, and the only relevant kinetic process is 
vibrational relaxation of the polyatomic molecules in the blend. 
For conditions characteristic of the use of substitute gases 
(i.e., moderate temperature and near atmospheric pressure), their 
vibrational relaxation time is on the order of 10 ”8 to 10” & 
seconds (Ref. 3), which corresponds to vibrational relaxation 
lengths very much smaller than crucial flow field dimensions 
along stream tubes. Equilibrium y values should therefore be 
valid for substitute gas flows. 

If the SCRJ combustor is efficient, the value of Dj- at the 
combustor exit should be of the order of unity because ti ^ 10~^s 
for nearly complete energy release in the supersonic diffusion 
flame (Ref. 4). As the exhaust flow expands, the local tempera- 
ture and pressure decrease , the reaction rate decreases (or reac- 
tion time increases), and Dj < 1. This case describes a condi- 
tion in which an essentially frozen composition experiences the 
expansion process. The detonation tube gases at the model test 
section can be tailored to provide the same composition as the 
flight SCRJ combustor exit condition, but at a much higher pres- 
sure to achieve proper scaling. The higher pressure results in 



more recombination reactions and tends to defer the onset of 
freezing. Calculations (Ref. 5) showed that this effect nearly 
compensates for the high dissociation in the stagnation region of 
the detonation tube. Although the chemical reaction history is 
not identical, the detonation tube system provides the closest 
reproduction of flight chemistry of any available simulation 
method. For the substitute gases, Dj » 1, and when used in 
SCRJ simulation facilities, they should always be in a thermody- 
namic equilibrium condition. 


The Dm parameter involves the relation of combustion 
energy release, and the ratio of residence to reaction time (or 
Dj) to the thermal energy transferred to the working fluid. Ex- 
amination of the test concepts we propose does not show any sig- 
nificant relevance of DjII to the SCRJ simulation process. 


The Eckert number is useful in characterizing the heat 
transferred from the SCRJ exhaust products to the afterbody sur- 
face and cowl, and thereby the effect of heat transfer on growth 
of the boundary layer. This parameter relates the prevailing 
temperatures of the exhaust products and the boundaries to the 
heat transferred to the walls, or 


Ec 



where 6 is the temperature difference between the exhaust gas 
at the outer edge of the boundary layer and the nozzle wall tem- 
perature. The SCRJ value for this parameter will not be dupli- 
cated directly in the detonation tube, unless the SCRJ wall is 
cooled to near room temperature, or the model is heated greatly. 
Fluid dynamic and fluid properties simulation is almost com- 
pletely reproduced, so Eckert number is reproduced only if wall 
temperatures are equal. We have concluded that the secondary 
effect of Ec on pressure distribution does not justify the 
large expense of heating the model wall and the resultant ef- 
fects on instrumentation. 


In substitute gas flow, the Eckert number can be matched to 
the SCRJ condition, despite the fact that d and the heat trans- 
fer is much smaller for the relatively cool test mixtures. This 
is possible because the model test gas velocity and Cp values 
are each quite different for the SCRJ gases at comparable model 
stations. Thus, for example, a model wall at room temperature, 
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using 400° K substitute gas, can match Ec for the SCRJ, where 
the exhaust gases are over 2000 ° K. 

It should be noted that the usual inviscid heat, transfer 
simulation, which involves only duplicating Nusselt number (Nu = 
qL(kA0)"-*- and geometry, is not adequate for the SCRJ exhaust 
nozzle process because of the large compressibility effects 
(i.e., large density changes) as well as viscous dissipation. 


Substitute Gases 


Substitute gases appear to be the best solution to the con- 
flict between the required thermodynamic behavior in afterbody 
model testing and the present lack of a method for generating high 
temperature exhaust gas chemistry and thermodynamics in a wind 
tunnel test. Low specific heat ratio, y, gases have been used 
with success for simulating high temperature air flows in hyper- 
sonic testing (e.g., Refs. 1, 6 , and 7), and those experiences 
were the initial stimulus to their application in SCRJ simulation 
Typical candidates for these substitute gases are several Freon 
compounds and sulfur hexafluoride (SF^) that can be used at 
moderate operating temperatures and pressures in a wind tunnel 
for supersonic speeds and are chemically stable over a wide tem- 
perature range. However, some of the best Freon compounds and 
SFfc require too low a temperature to duplicate the specific heat 
ratio of SCRJ gases at the entrance of the exhaust nozzle, and, 
consequently, they experience phase change (condensation) along 
the expansion process. This undesirable test gas property can be 
overcome by fractional additives of high specific heat ratio, 
stable noble gases, such as argon. We thereby create binary and 
ternary mixtures capable of matching the ratios of thermodynamic 
properties of virtually any SCRJ high speed flight propulsion sys 
tem (exemplified by the values of y of the combustion products) 
Figure 1 shows the y(T) behavior of several such blends as well 
as that of the SCRJ exhaust flow in equilibrium. It should be 
noted that we do not wish to match y(T) ; instead, we wish to re 
produce y(X) at a much lower temperature level, which can be 
viewed in the previous context as matching y (T/T 3 ) . 

Thus, by selection of appropriate and convenient working 
temperature, pressure, and binary gas mixture composition, the 
SCRJ exhaust distributions of y, M, and Re can be matched in 
a test facility using small scale models. Furthermore, the dif- 
ferent working requirements of different test techniques can be 


Fig. 1 Variation of y 
SCRJ Exhaust i 
Formulations 






accommodated. Figure 2 presents curves of how changes in compo- 
sition affect the temperature at which there is an exact match 
to y and M at the exit for flight conditions of = 8, 6, 
and 4; an effective angle of attack (a + 3) of 4°; and 
equivalence ratios (<t>) of 1.0, 1.0, and 0.6, respectively. 
These representative trends are for binary mixtures of Freon 12 , 
13, 13B1, and 14 with argon additive and can readily be extended 
to match requirements of other flight conditions or for other 
matching points in the flow. 

Table 1 compares several similarity parameters and thermo- 
dynamic and transport properties of a M^ = 8 SCRJ operational 
condition with two typical substitute gas mixtures of Freon 12 
and 1351 with argon, additive . Figure 3 shows that these mix- 
tures are designed to produce the best pressure distribution 
match to the 1-D computed SCRJ exhaust gas along the nozzle 
walls. Small compositional changes in the selected substitute 
gas mixture can result in significantly lower stagnation tem- 
perature requirements, with relatively minor alteration in pres- 
sure distribution match. Other flight Mach numbers can also be 
matched similarly. It is necessary to prevent the Freon com- 
ponent gas from condensing within the model flow field. To il- 
lustrate this point, the required 622 K stagnation temperature 
for the 50 percent Freon 13B1 + 50 percent Argon mixture (see 
Table 1) can be lowered to 485 K by changing to a 58 percent 
Freon 13B1 & 42 percent Argon mixture. The Freon fraction then 
will still be well within the vapor phase — a static tempera- 
ture of 185 K and partial pressure of 4500 Newton/m^ (0.653 
psia) — — for Mach 4.06 flow at the end of the model afterbody. 
The 22 percent lower stagnation temperature can be extremely 
effective in simplifying the utilization of existing wind tun- 
nels. Computational methods are described in the appendix. 
Subsequent calculations discussed in the Two Dimensional After- 
body Pressure and Force Predictions section will show how two- 
dimensional calculations using the same substitute gas behavior 
compare to the SCRJ case . 

The SCRJ exhaust gas properties and similarity parameters 
are obtained from the standard NASA Lewis Research Center Com- 
puter Program (Ref. 8) for an equilibrium kinetics expansion 
process of hydrogen/air reaction products. The program results 
are augmented by our additional calculations of Sc and Le , 
which use a mutual diffusivity value computed for the SCRJ ex- 
haust gases. This value was obtained by assuming that the 
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Fig. 3 Preliminary Match of SCRJ Exhaust Pressure Distribution 
with Five Substitute Gas Expansions 



principal SCRJ exhaust components of N 2 and H 2 O can be ap- 
proximated as a binary mixture. Calculations of the substitute 
gas properties and parameters are described in the appendix. 

The properties of the two substitute gas flows shown in 
Table 1 were established by matching the Mach number and Rey- 
nolds number exactly for the conditions at the combustor exit 
of a f the scale model once the gas blends were chosen. The 
gas flows were all assumed to be homogeneous. The initial 7 
of each of the substitute gases is larger than that of the SCRJ, 
but these selected binary mixtures produce nearly the same over- 
all pressure distribution as the SCRJ along the nozzle expan- 
sion. As was shown in Fig. 1, the y versus T trends of the 
SCRJ gas and substitute test gases are not identical; however, 
it is possible for the latter to reproduce the working range of 
y for the SCRJ, but at much lower temperatures. As shown by 
Fig. 3, the matching of the different substitute gas isentropes 
to the SCRJ expansion pressure distributions is very close. 

In preliminary selection of substitute gas candidates, the 
best similarity parameter to be used is the value of y charac- 
terizing the over-all expansion process. For isentropic expan- 
sion of the chemically nonreactive substitute gas; mixtures, con- 
sidered as a perfect gas, the process-averaged y is the isen- 
tropic exponent of the relation (P/p)^ = constant. Then 

_ log (P 2 /P 1 ) log (P 2 /P 1 ) 

7 ” log (pJTp'JT " log (P 2 / p i) " lo s (t 2 / t i) 

where the subscripts 1 and 2 refer to the nozzle entrance 
and exit conditions, respectively. For the three candidate sub- 
stitute gas compositions whose gas dynamic properties are given 
by Table A-l (see p. 68 ), the 7 for the nozzle expansion pro- 
cess (between Mach numbers 2.36 to 4.06) is as follows: 

50 percent Freon 13B1 + 50 percent Argon - 7 = 1.24 

40 percent Freon 12 + 60 percent Argon - 7 = 1.27 

58 percent Freon 13B1 + 42 percent Argon - 7 = 1.20 

The expansion of SCRJ gases involves small changes in equilib- 
rium chemical composition leading to slight changes in molecular 
weight dependent on the local pressure and temperature. In this 
polytropic process, the process-averaged exponent, 7 , governing 
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-»c 

the (P/p) y relation is now defined as 

-c L °g <yy 

7 log (Pj/Pj^) - log (T 2 «! 1 /t 1 W 2 ) 

to account for the molecular weight changes. From the equilib- 
rium gas dynamic data generated by the method of Ref. 8 for the 
M m = 8, a = 4 , cp = 1.0 SCRJ operational condition, the 

value of y equals 1.24. For the frozen composition SCRJ gas 

dynamic properties the value of y is calculated by the same 
relation as for the substitute gases, and for the M^, = 8 SCRJ 
condition has a value of 1.277. From data shown in Fig. 16 it 
will be seen that the resulting forces produced by the SCRJ 
pressure distribution along the nozzle surface is nearly exactly 
simulated fry substitute gas expansions where y is very nearly 
equal to y , in the SCRJ equilibrium case, and where the y's 
are equal for the frozen SCRJ case. Note again that the purpose 
of y is only for preliminary selection; the actual y(t/t 3 ) 
function for the substitute gas must be used in a multidimen- 
sional calculation before it can be accepted as a valid substi- 
tute gas. 

Other similarity trends shown by the data of Table 1 are: 

• Prandtl No. The low value (Pr ~ 0.50) for the 
SCRJ exhaust shows that thermal diffusion is stronger 
than momentum diffusion. The substitute gas selec- 
tions will ha\e a somewhat less dominant thermal dif- 
fusion process, being more like the familiar levels 
of low temperature air. In this regard, the Freon 12 
and argon mixture is a better simulation than the 
Freon 13Bl-argon composition. Note that Pr is not 
changed significantly by pressure changes, nor does a 
change in Pr have a large effect on pressure dis- 
tribution in nearly adiabatic attached flows. 

* Schmidt No. In view of the large Sc value for 
both the SCRJ and the substitute gases, the mass dif- 
fusion in each case is £ to .■§ of the diffusion of 
momentum in the flow field. All cases are comparable 
although the Freon 12 + argon component mixture is a 
slightly better match to the SCRJ case. 
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Lewis No. 


Lewis No. The very low value of the SCRJ and sub- 

stitute gases indicates that in both situations, the 
diffusion of the temperature field is almost an order 
of magnitude greater than mass diffusion. 


As previously discussed, dynamic similarity requires simu- 
lation of gas properties; for detailed nozzle boundary layer de- 
velopment the simulation of transport property variation with 
temperature is of great importance because of the large tempera- 
ture gradients in the gas adjacent to the nozzle wall. While 
the Reynolds number matching of model and SCRJ governs the over- 
all flow process, the boundary layer bears special consideration 
because of its possibly significant interaction with the flow 
field pressure distribution. The viscosities of the binary sub- 
stitute gas mixtures listed in Table 1 are about one third that 
of the SCRJ exhaust gas at the nozzle entrance. This reflects 
primarily the large difference in temperature of the two working 
fluids . 


Perhaps of greater significance is the temperature varia- 
tion of viscosity of the substitute gases compared to the SCRJ 
gases in their respective temperature ranges of operation as 
shown by Fig. 4. . This variation can be characterized by the 
value of the exponent, n, in a power law approximation of p, 
versus T. As shown in Table 1, the n value for the binary 
substitute gases is about 10-12 percent larger than the SCRJ 
gases. 

It is well-known that the viscosity-temperature exponent 
will have a greater effect on boundary layer growth and related 
parameters at supersonic Mach number than will relatively large 
changes in Pr. Because the inviscid similitude laws provide a 
correspondence in temperature ratios between corresponding adia- 
batic points in substitute and SCRJ flows [i.e., (T 2 /T 1 ) sub- 

stitute gas ^ (T 2 /T]) SCRJ], we can establish reference vis- 
cosity ratios for the afterbody boundary layer characterized by 
adiabatic temperature levels — say, stream and stagnation tem- 
perature. These will ensure corresponding viscosity distribu- 
tions in the wall cooling (Eckert number) are all the same. It 
is our feeling that the matches of n and Pr shown in Table 1 
are more than adequate for their relatively weak coupling to 
pressure distribution, but the final test of this opinion will 
come when substitute gases are tested against detonation tube 
SCRJ flows in Phase II of this program. 



LEGEND : 





The viscosity properties of a ternary mixture are also 
shown by Fig. 4 where a Freon 13Bl-sulfur hexafluoride -neon 
mixture formulation is plotted. 

The ternary mixture was examined to see if its y(T) vari- 
ation would prove superior to the best binary mixtures, or if 
there were significant advantages in the transport properties of 
the ternary mixture. The limited data we generated on ternary 
mixtures did not indicate any significant advantage for them; in 
fact the (X - T exponent moved further from that of the SCRJ 
exhaust gas. 

A slightly better (i.e., lower) value of n could be 
achieved by using a substitute gas mixture of Ne and Freon 
13B1 at a stagnation temperature of about 1600 K, While this 
could be easily done in a shock tunnel, it would create opera- 
tional difficulty and expense in the subsequent wind tunnel 
testing that would preclude its selection. 

In summary, the distribution of transport property effects 
in the SCRJ exhaust and binary mixture substitute test gases 
follow the same relative order of dominance, namely 

thermal > momentum > mass , 

and very similar viscous aerodynamic behavior should be encoun- 
tered within each working medium. 


Mixing Layer Similitude 

Complete SCRJ exhaust simulation eventually will involve 
the mixing of simulated exhaust and external flow fields. The 
scaling laws for turbulent mixing are not known reliably, so 
either all of the suspected important scaling parameters must be 
matched or an experimental investigation must be designed in 
which any imperfect matching can be varied and evaluated for its 
effect. The major mixing similitude parameters are: Re, ini- 

tial boundary layer properties on the internal and external sur- 
faces (perhaps scaled to the inviscid flow properties), and the 
ratios of internal to external density, velocity, viscosity, and 
temperature. For inviscid similitude, a simple analysis indi- 
cates that the parameter (T^/T e ) Ql\ e /^j) > where 771 is molecu- 
lar weight, will simulate some of these mixing properties if 
maintained equal in the model and prototype. (The subscripts i 
and e refer to internal and external flows, respectively.) 


For example, the density ratio across the mixing layer is: 




,p. 


n.. 


and the velocity ratio is: 

V 


M.. . 

V \M A-v T IK.) 


e e 1 


For the SCRJ, the magnitude of the parameter (T-p/T e )( at 
the cowl lip is approximately 8.5. This value should be matched 
in the detonation tube simulator tests and in wind tunnel tests 
where relatively low temperature substitute gases are used. 


For the detonation tube, the ratio (T-p/^-p) is about 80 
so the required (T e /% e ) ratio is about 9.5 for exact match- 

ing. This condition is easily obtained by using air (7K e = 29) 
for the external flow, and setting test conditions to produce a 
T e of about 275 K at the cowl lip station. 


For the wind tunnel test situation, the types of substitute 
gases selected for SCRJ internal flow simulation give a (TjV^i) 
ratio of 3.7 to 4.7. This establishes a (T e //ft e ) ratio re- 
quirement of 0.43 to 0.55 for the external flow substitute 
gas. For the safe Freon compounds and binary mixtures considered 
for SCRJ exhaust gas simulation, the molecular weight will be in 
the range of 70 to 150. If used as an external gas, this 
would result in a need for extremely low temperature, T e (on 
the order of 50 K to 100 KX at the cowl lip location with its 
accompanying danger of phase change for the Freon component at 
practical pressure levels. To obtain more practical temperature 
and pressure levels a different class of fluorocarbons is needed. 
A preliminary survey of safe, chemically stable candidates for 
this purpose indicates that the Freon "E" series of high molecular 
weight, low 7 , compounds, with (T e /^! e ) values of 0.52 to 
0.72 at their boiling points could be useful components of a bi- 
nary mixture with high molecular weight noble gases. Freon E5, 
for example, has a molecular weight of 950 and a boiling point, 
at atmospheric pressure, of 497 K. Using these "E" series 
fluorinated ethers at lower pressures than atmospheric could re- 
duce the temperature requirements and decrease the (T e A e ) ratio 
to the full extent needed to examine mixing processes in wind tun- 
nels. Clearly, this problem requires more extensive study which 
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could be pursued in the context of a research program concerned 
with mixing phenomena of high speed parallel flows of dissimilar 
gases. Also requiring further consideration are the methods and 
costs of implementation by which the external and internal flows 
can be supplied and sequenced in short duration tests in both 
the wind tunnel and detonation tube simulator test techniques. 
However, a conceptual solution for substitute gas simulation of 
SCRJ mixing processes appears available within the framework of 
present technology. 
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DETONATION TUBE SIMULATOR 


The Grumman detonation tube simulator to be used for the mea- 
surement of exhaust flows and validation of the substitute gases 
is capable of giving nearly correct chemistry and total enthalpy 
for the hydrogen/air combustion system throughout the entire pro- 
posed flight regime. Reynolds numbers can easily be maintained 
at flight values by running at elevated pressures with both com- 
bustion and substitute gases. Comparisons of pressure distribu- 
tions over a model afterbody obtained with the combustion gases 
against those obtained with candidate substitute gases will be 
used to select the proper substitute gas and operating conditions 
for subsequent wind tunnel tests so that the necessary dynamic 
similitude prevails. 

A schematic of the detonation tube facility is shown in 
Fig. 5. It consists of a 6.1 meter long, 7.62 cm inside 
diameter driver section initially separated from a 10 meter 
long, 12.7 cm inside diameter driven tube by a metallic dia- 
phragm. The driven tube is terminated by a supersonic nozzle de- 
signed to produce a flow that will match the expected exit plane 
conditions of the proposed scramjet engine. The nozzle exhausts 
into a 1.83 meter diameter, 3.66 meter long test section. 

Mounted within the section, and mating to the supersonic nozzle 
is the model afterbody. For the proof -of-concept experiments, 
the upper afterbody surface will be a flat plate instrumented 
with up to 20 pressure transducers. The test section can be 
evacuated prior to a run to any desired pressure down to 
5 x 10 ’5 torr. The facility has been designed to handle com- 
bustible hydrogen mixtures safely. Various interlocks and leak 
detection devices are incorporated into the automatic gas handling 
system. Details of this system may be found in Ref. 1. 

The facility has available up to 35 channels to carry model 
data. All data are recorded through FM/FM multiplexing on mag- 
netic tape and reduced by a special off-line computer. Heat 
transfer rates are measured using thin film platinum resistance 
gauges. Pressures in the past have been measured with standard 
piezoelectric gauges for pressures above 7000 N/m 2 or with spe- 
cial Grumman designed gauges for low pressure work. For the 
afterbody pressure measurements above 7000 N/m 2 we intend to 
use a new system using unbonded semiconductor strain gauges. 

Those are much lower impedance devices than the piezoelectric 
type and will lead to much less mechanical noise. These instru- 
ments use a bridge circuit; differential amplifiers can therefore 




Fig. 5 Detonation Tube Simulator Schematic 




be used which will allow the elimination of common-mode noise. 
Other facility induced noise problems will be unknown until we 
actually try the system. Since capability for accurate pressure 
measurement is critical to the successful conclusion of the next 
phase of the program, we are devoting considerable effort outside 
the coverage of this contract to achieving the maximum instrument 
capability possible. For pressures above 7000 N/m^. gauges are 
available with sufficient frequency response and with accuracy to 
- 5 - of one percent of full scale. In the pressure range from 
70 N/m2 to 7000 N/m^ we are currently investigating the use of 
microphone- type gauges to obtain similar accuracies. 


Photographic coverage available includes time-exposure glow 
photographs in color or black and white, high speed movies, and, 
where there is sufficient density, schlieren photographs. We hope 
that the schlieren photographs will be sufficiently sensitive to 
measure relative boundary layer growth in substitute gas and de- 
tonation tube tests. 

Two modes of detonation tube operation are possible. In the 
first, called the incident or forward-running detonation tech- 
nique (Ref. 9), the driver section is pressurized with helium to 
a pressure high enough to insure detonation of the combustible 
gases when the main diaphragm ruptures. The diaphragm is scribed 
to a predetermined depth so that it will rupture at a specific 
pressure. The shock wave formed will rapidly become a self- 
propagating detonation wave that will reflect from the nozzle end 
of the driven tube as a shock wave. The gas is heated and pres- 
surized by the combustion behind the detonation wave and further 
by the energy addition behind the reflected shock wave. A weak 
diaphragm at the entrance to the nozzle is ruptured by the arrival 
of the incident wave and the high energy slug of gas expands 
through the nozzle into the test section. In order to prevent the 
rapid expansion that exists behind a free-running detonation wave, 
the initial helium driver gas pressure is chosen so that when it 
expands after the rupture of the diaphragm its pressure matches 
that immediately behind the detonation wave. This produces quasi- 
steady flow conditions for two to four milliseconds* 

The second mode of operation is called the backward-running 
detonation technique. This technique was developed specifically 
for the simulation of hydrogen/oxygen rocket engine plumes (Ref. 1) 
where the incident technique would have produced an enthalpy higher 
than the actual engine case because of the shock wave heating added 
to the heat of combustion. A very weak shock wave is sent into the 
combustible mixture in this method (through the use of a low driver 
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pressure and weak main diaphragm) . The idea is to prevent igni- 
tion on the incident wave, but have the gas ignite when the in- 
cident wave reaches the nozzle end of the tube. The detonation 
wave then travels away from the nozzle back towards the driver 
section. The expansion following this wave drops the pressure 
and enthalpy in the test gas, which then expands through the 
nozzle. We have found in previous work (Ref, 1) that the igni- 
tion at the end of the tube is best accomplished by a timed 
electrical spark, although under certain conditions spontaneous 
ignition behind the reflected wave will occur. 

The choice of technique to be used for the scramjet simula- 
tion depends on the final conditions required and to some extent 
on operational experience during exploratory tests to be made 
during the next phase of this program. Further discussion of 
this choice is given in the next section where calculated running 
conditions for various mixtures are presented. 

For running the inert substitute gases in the detonation 
tube we run the facility as a conventional reflected shock tun- 
nel. That is, the test gas mixture is placed in the driven tube. 
Helium driver gas is used. The shock wave formed on rupture of 
the main diaphragm reflects from the nozzle end of the driven 
tube to produce a slug of high energy gas which expands through 
the nozzle into the test section. The stagnation conditions of 
the substitute gas are determined by the strength of the incident 
shock wave in the driven tube, which may be easily controlled. 


Detonation Tube Simulator Running Conditions 

We have been able to show that by the use of one of our two 
detonation tube techniques we will be able to simulate the HRA 
exhaust nozzle flow with exact total enthalpy (flight and com- 
bustion) , full scale Reynolds number, and nearly exact chemistry 
(exact elemental composition, with nearly exact molecular compo- 
sition) . Calculations were made for the three basic flight con- 
ditions specified in NASA Langley Statement of Work 1-53-3387, 
March 12, 1973, plus the M TO = 8 flight conditions. It appears 
at this stage that the detonation tube techniques will give ex- 
act total enthalpy simulation only up to Moo = 10 at i> = 1, 
but close simulation is available to higher flight Mach numbers. 
In the following pages we give some details of the calculations 
which show the broad range of simulation capability of the deto- 
nation tube and how the above mentioned flight conditions fall 
within our operating range. 
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Total Enthalpy and Chemical Composition 

With the backward-running detonation tube technique we can 
give exact total enthalpy and equilibrium chemical composition 
simulation for M m = 4 at <t> = 0 . 6 , = 6 at <t> = 1 , and we 

may use this method as high as M*, = 8 at <J> = 1 . 


The procedure for determining the tunnel stagnation condi- 
tions is as follows (see also Ref. 1). For a given incident wave 
speed at the end of the tube, the state properties and particle 
velocity behind the wave are calculated, assuming real-gas 
thermodynamic properties and frozen composition. The reflected 
wave is assumed to be a Chapman -Jouguet (C-J) detonation whose 
initial conditions are the conditions behind the weak incident 
shock wave. The unsteady expansion from the C-J state to the 
tunnel stagnation conditions is calculated assuming an isentropic 
expansion with equilibrium thermodynamic properties. The govern- 
ing equation for this expansion is a - U = constant, where 


P P 


a 


(a/p) g dp . 

p ref 


For the C-J state, U = (Ucj - U 2 ) > and for the end wall 
state (Region 5), U s 0. The chemical equilibrium program of 
Ref. 8 is used throughout these calculations. 


There are a wide variety of gas mixtures and initial shock 
wave speeds that we could use to achieve these simulations. Mix- 
tures would consist basically of H 2 , 02 > N 2 , and Ar, with vari- 
ous amounts of NH 3 or N 2 O added to achieve the desired en- 
thalpy. For 4> = 0.6 the gas mixtures could be (Y)NH 3 + 

(0.2511 - 3/2 Y)H 2 + (0.2095)02 + (0.7809 - 1/2 Y)N 2 + (0.0096)Ar 
with 0 < Y < 0.1674, or they could be (X)N 20 + (0.2511)H2 + 
(0.2095 : 1/2 X)0 2 + (0.7809 - X)N 2 + (0.0096)Ar with 
0 < X < 0.4190. 

For <t> = 1 the mixtures could be (Y)NH 3 + (0.4190 - 
3/2 Y)H 2 + (0.2095)0 2 + (0.7809 - 1/2 Y)N 2 + (0.0096)Ar with 
0 < Y < 0.27933, or (X)N 2 0 + (0.4190)H 2 + (0.2095 - 1/2 X)0 2 + 
(0T7809 - X)N 2 + (0.0096)Ar with 0 < X < 0.4190. 

The actual gas mixture used will be chosen for operational 
convenience in obtaining the desired stagnation pressure, P 5 . 



Figures 6 and 7 show the variation of total enthalpy with initial 
shock wave speed for these gas mixtures. The problem of achiev- 
ing the Moo = 8 flight case with the backward-running technique 
occurs because the incident shock wave Mach number must be greater 
than 2.5 (see Fig. 7 ). If this high speed causes an incident 
detonation, we will use the front-running technique to simulate 
M^ = 8 . The backward-running technique is preferable because it 
requires a much lower driver gas pressure to achieve the same 
stagnation pressure. The higher driver pressure required by the 
front-running technique is due to the high contact surface ve- 
locity required to match the pressure behind the detonation wave. 

With the front-running detonation tube technique we will be 
able to achieve exact total enthalpy simulation at $ = 1 from 
about Moo = 7 to M^ = 10, overlapping the conditions simu- 
lated by the backward -running technique. For a given gas mixture 
the forward-running technique generates a higher total enthalpy 
because the burned gas is further heated by a reflected shock 
wave, rather than cooled by an expansion wave as with the back- 
ward-running technique . 

The stagnation conditions for each gas mixture are determined 
by first calculating the conditions behind the detonation wave 
assuming the C-J condition (see Ref. 9 ). The steady flow con- 
servation equations and the equation of state are then applied 
across the reflected shock wave assuming a constant ratio of spe- 
cific heat, y, equal to that behind the detonation wave. 

Figure 8 is a graphical presentation of the variation of 
total enthalpy with the <t> = 1 gas mixtures specified at the 
beginning of this subsection. For a specified total enthalpy, 

X and Y are only weak functions of the stagnation pressure, 

P5, since the initial shock Mach number is now fixed at the pre- 
determined Chapman- Jouguet detonation wave speed for each mixture. 
It is fortuitous that to simulate M ro = 8, the flight condition 
agreed upon for our initial experimental work, no other molecular 
forms of the elements in the combustion gases are required to be 
added to the basic H2/air mixture (see Fig. 8). 


Model Size and Reynolds Number 

Both the backward -running and forward-running techniques re- 
quire reflected shock waves. As a general rule, to produce the 
desired stagnation conditions behind a reflected shock wave, the 
throat area in the end wall of the shock tube should not exceed 
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10 percent of the shock tube cross section area. The combustor 
exit plane Mach numbers and shock tube nozzle exit to throat 
area ratios required to generate those Mach numbers are given 
below for the cases under consideration. 


Shock Tube 


M 

00 

4> 

M 

exit 

A e /A* 

4 

0.6 

' 1.29 

1.07 

6 

1.0 

1.68 

1.40 

8 

1.0 

2.34 

3.10 

10 

1.0 

3.00 

8.30 


Note: The variation of Me X it anc * A e / A with a + P is 

very small. The above numbers are average values 
for the 3 a + P's at each Moo- 

If we assume a two dimensional combustor with a square exit 
plane and a throat area A* equal to one-tenth that of the de- 
tonation tube cross section, the maximum exit plane dimensions 
become : 


M Maximum Exit Plane 

°° Dimensions, cm 

4 3.683 x 3.683 

6 4.318 x 4.318 

8 6.1 x 6 . 1 

10 10.16 x 10.16 


A rectangular exit plane with an aspect ratio greater than 
1 was chosen, primarily because it will give a relatively larger 
volume of two dimensional flow than the alternatives, giving 
measurements which could be more easily compared to theory. The 
M^ = 8 flight condition was also chosen as the nominal baseline 
for the first series of experiments required to prove the con- 
cept of detonation tube simulation for a SCRJ. In Fig. 9 we show 
the maximum combustor exit plane dimensions for various aspect 
ratios for the four flight conditions under consideration. Given 
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the restraints of a 12.7 cm i.d. detonation tube with throat 
area one tenth of the detonation tube area and = 8 , we chose 
a 3:1 aspect ratio for the combustor exit plane with a width 
of 11.43 cm and a height of 3.81 cm. With the 3.81 cm 
height we can also fit three 3.81 cm x 3.81 cm square com- 

bustor nozzles in the end wall of the detonation tube, thus simu- 
lating the 3:1 aspect ratio with and without separators. Fig- 
ure 10 is a sketch showing the proposed single nozzle arrangement 
and Fig. 11 shows the same shock tube and afterbody model, but 
with the three smaller nozzles installed in the detonation tube 
end wall. 

To achieve full scale Reynolds number we first determined 
the exhaust flow Reynolds number per meter at the prototype com- 
bustor exit plane, using the NASA/Lewis computer program (Ref. 8 ). 
These results are presented in Table 2. Since we chose a model 
combustor height of 3.81 cm while the prototype is assumed to 
be 30.48 cm, we will be doing 1/ 8 scale testing based on 
combustor exit height and must achieve in the model combustor a 
Reynolds number per meter eight times that of the prototype for 
proper viscous flow development on the afterbody. 

To determine the stagnation conditions in the detonation 
tube that are required to achieve the desired Reynolds number at 
the model combustor exit, we assumed the following: 

1. The total enthalpies of the prototype and model 
flows are identical. 

2. The equivalence ratios of the prototype and model 
flows are identical. 

3. The flow from the detonation tube stagnation region 
to the model combustor exit is in equilibrium. 

4. The combustor exit to sonic throat area ratio in 
the model (A 3 /A*) is the same as the theoretical 
combustor exit to sonic throat area ratio in the 
prototype . 

The last assumption required calculating a theoretical sonic con- 
dition for the prototype. The differences in exit Mach number, 
between model and prototype, were on the order of a few hundredths 
of a Mach number, not enough to justify additional calculations to 
achieve an exact match. Using these assumptions with the computer 
program of Ref. 8 , we determined the detonation tube stagnation 



Detonation Tube 



Fig. 10 Sketch of Detonation Tube and Model Installation 

for Unsegmented, = 3 Scramjet Combustor Exit Plane 




EQUILIBRIUM THERMODYNAMIC AND TRANSPORT PROPERTIES AT SCRJ COMBUSTOR EXIT PLANE 
WITH CORRESPONDING ISENTROPIC STAGNATION CONDITIONS 
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Note; 1 atm - 1.01325 x 10 ° N/tn 



required to achieve the desired Reynolds number, together-with 
the other thermodynamic and transport properties at the model 
combustor exit plane. 

Table 3 summarizes the results of these calculations for the 
Mo,, = 8 flight cases. The differences between the 1/8 scale 
model and prototype exit plane Mach numbers and other state 
properties may be seen by comparing Tables 2 and 3. Achieving a 
stagnation pressure between 60 and 110 atmospheres with back' 
ward-running detonations is no problem in our facility. We have 
repeatedly run different simulated propellants at pressures be- 
tween 10 and 400 atmospheres. The driver pressure, P 4 , re- 
quirements for these backward-running detonations are not at all 
severe — 10 to 20 atmospheres at the most. 

To achieve stagnation pressures, P 5 , of several hundred 
atmospheres with front-running detonations, higher driver pres- 
sures are required. For the Mach 8 flight cases using room tem- 
perature helium as the driver gas, the driver pressures will be 
between 200 and 300 atmospheres. We have routinely run for- 
ward-running detonations with P^/s between 300 and 400 at- 
mospheres and have the capability to operate at double that 
level. Thus, simulating the M^ = 8 flight case with 1/8 scale 
models at full scale Reynolds number is well within the operating 
range of the Grumman Detonation Tube Simulator. 

Ensuring that the edges of the exhaust expansion regions 
will be properly configured requires that the contact surfaces 
between external and internal gas in the detonation tube test lie 
outboard of the contact surfaces in the flight and wind tunnel 
cases. This is achieved by ensuring that the test chamber pres- 
sure in the detonation tube is well below the static pressure of 
the flight contact surface (i.e., the pressure at which internal 
and external flows can turn to the same flow direction) . The 
domain of validity for the test is then determined analytically 
for each flight condition. 


Experimental Verification of Simulation Technique 

The next phase of this program will include proof -of-concept 
experiments to validate the proposed technique. During these 
experiments measurements of the wave speeds and stagnation pres- 
sure, P 5 , will be made to insure that the desired stagnation 
conditions have been achieved. Calibration of the flow from the 
combustor exit will be accomplished by stagnation pressure and 
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TABLE 3 


EQUILIBRIUM THERMODYNAMIC AND TRANSPORT PROPERTIES AT f-SCALE 
MODEL COMBUSTOR EXIT PLANE FOR SIMULATING FULL SCALE REYNOLDS 
NUMBER WITH CORRESPONDING DETONATION TUBE ISENTROPIC STAGNATION 

CONDITIONS, M =8.0, 4> = 1.0 

00 \ 


a 

f p (deg) 

4 

8 

12 

P 3 

(atm) 

4.2731 

6.1788 

8.0730 


(°K) 

2310.0 

2316.0 

2368.0 

V 3 

(m/sec) 

2317.28 

2316.93 

2272.87 

M 3 

* 

2.386 

2.380 

2 .311 

V A 


3.1233 

3.0833 

2.8142 

Rey . 

No. ^ (ni ) 

1.736 x 10 7 

2.500 x 10 7 

3.081 x 10 7 

7 3 


1.2045 

1.2077 

1.2051 

’S 


24.524 

24.536 

24.522 

^3 

(p. poise) 

738.0 

739.0 

752.0 

ko 

p. cal n 
(cm-sec-°K/ 

f cal n* 

642.0 

622.0 

650.0 

3 

c 

0.5183 

0.5073 

0.5167 

P 3 

pr 3 

\gm-°K/ 

0.5953 

0.6029 

0.5979 


DETONATION 

TUBE ISENTROPIC 

STAGNATION CONDITIONS 

P 5 

(atm) 

66.96 

95.623 

109.45 

T 5 

(°K) 

3365.0 

3397.0 

3409.0 

H q 


578.4 

578.4 

578.4 

5 

a 5 

^gm / 
(m/ sec) 

1173.3 

1178.7 

1180.8 

^5 


1.1636 

1.1665 

1.1676 

*5 


23.651 

23.712 

23.735 


Note: 1 atm = 1.01325 x 10^ N/m^ 



heat transfer surveys using instrumented rakes at several sta- 
tions along the afterbody. A limited amount of spectroscopic/ 
radiometric analysis for exhaust composition and temperature 
assessment will also be conducted. 

Calculation procedures proven in previous plume simulation 
work (such as reported on in Ref. 9), will be used for compari- 
son with the measured data. We have shown in the past that such 
comparisons give a valid indication of the fidelity of a given 
detonation tube simulation. In addition, IR measurements ob- 
tained in plumes simulated by the detonation technique have been 
shown to give remarkable agreement with other experimental tech- 
niques as well as with detailed theoretical calculations (Ref. 
10). The methods used in such theoretical computations are pre- 
sently limited to axisymmetric nozzles, so their application to 
the nonuniform three dimensional afterbody flow of this project 
will only be an intermediate step. Final and complete theoreti- 
cal substantiation of our data in the flight situation must 
await the completion of theoretical methods now under develop- 
ment by other contractors. 


TWO DIMENSIONAL AFTERBODY PRESSURE 
AND FORCE PREDICTIONS 


The goal of this program is to provide a system of substitute 
gases for use in a wind tunnel model that will duplicate the dis- 
tribution of forces produced by the scramjet engine exhaust over 
the vehicle afterbody. In this section we present calculations 
of pressure distributions for several different cases, and, in 
addition, show integrated axial and normal force values. All cal- 
culations in this section are for a flight Mach number of 8. 

The results are the basis of our conclusion that it is possible 
to select a substitute gas to fulfill the goal stated above. 

We have performed two dimensional method of characteristics 
calculations for nine separate cases. For the full scale flight 
hydrogen/air combustion system, we did four calculations: equi- 

librium and exit plane frozen, each at inlet ramp angles (a + p) 
of 4° and 12° . We have similarly calculated pressure distri- 
butions for an eighth scale detonation tube simulation, matching 
all conditions at the exit plane to that of full scale, but 
raising the pressure to duplicate Reynolds number. The eighth 
scale calculations were carried out for a . + p of 4° and 12°, 
but only for equilibrium chemistry. An exact nonequilibrium two 
dimensional kinetic calculation has not been carried out. 


Three substitute gases were chosen using the methods de- 
scribed earlier to match the full scale case of a + p = 4° and 
equilibrium chemistry. The first consisted of 50 percent by 
volume Freon 13B1 and 50 percent argon, the second 40 percent 
Freon 12 and 60 percent argon, and the third 58 percent Freon 
13B1 and 42 percent argon. Mach number and Reynolds number at 
the nozzle exit were matched for the substitute gas calculations. 
We also calculated a pressure distribution assuming y = 1.4 air 
with the same exit Mach number as an illustration of the sensi- 
tivity of loads to large errors in gas characteristics. 

In all cases the two dimensional method of characteristics 
computer program of Ref. 11 was employed. The thermodynamic 
variables for the combustion gases were generated by the computer 
program of Ref. 8, while those for the substitute gases were cal- 
culated as described in the appendix. The flow properties at the 
combustor exit were assumed uniform, parallel, and in equilibrium 
for all cases. 


48 



The geometries of the afterbody nozzle surface and the ccwl 
surface are described by the following equations: 

Afterbody 

Y » 0.7187 X 2 + 1.0 0 < X < 0.25 

Y = 0.35935 X + 0.95508 0.25 < X < 18.54 

Cowl 

Y = 0 0 < X < 1.11 

Y » -0.4204 X 2 + 0.933 X - 0.518 1.11 < X < 1.235 

Y = -0.1051 X + 0.1232 1.235 <X < 3.12 

where Y = Y/Yo, X = X/Yo, and Y 3 is the height of the nozzle 
exit. 

A sketch of this geometry showing the origin of the coordi- 
nates is shown in Fig. 12. Also shown in this figure is the 
characteristic diagram calculated for the full scale engine, for 
M^ = 8 , a + p = 12°, and assuming fully equilibrium flow. The 
flow fields for all other cases mentioned above with the exception 
of the air case differed from the one shown only in small detail. 
Regions of uniform flow are labeled "U." At the cowl lip there 
is a Prandtl-Meyer expansion to match the local external pressure. 
In all of the cases considered except the 7 = 1.4 case the first 
characteristic line from the lip intersected the upper surface at 
X ~ 15.5. The characteristic line that intersects the downstream 
end of the upper surface corresponds to an expansion of about 3° 
(through a pressure ratio of 1.15). It has been shown that the 
turning pressure ratio at the cowl lip is close to 2 for the 
prototype flight case. This point is very important, as it in- 
dicates that at least for two dimensional testing the internal 
flow is all that need be simulated to ensure that the afterbody 
ramp pressures are valid. Mixing effects are not likely to de- 
stroy this independence of the external flow as long as the com- 
bustor exit is everywhere underexpanded, but three dimensional ef- 
fects at the ends of the combustor will require more complete 
simulation for regions near and outside the diverging lateral 
contact surfaces. 



Fig. 12 Nozzle Afterbody Geometry and Characteristic Diagram 
for Full-Scale Engine at M«>= 8 ,a+0 - 12° 
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In Fig. 13 we have plotted the pressure distributions, non- 
dimens ionalized by the combustor exit plane pressure for each 
case, for the combustion gases at exit conditions corresponding 
to a + p = 4° . In Fig. 14, we show pressure distributions for 
a + P = 12°. It. can be seen that the pressures calculated as- 
suming fully equilibrium flow are slightly higher than those cal- 
culated using frozen chemistry. Also the excellent agreement 
with the case calculated for the detonation tube simulation is 
evident. In Fig. 15 we have compared the calculated pressures 
for the eighth scale detonation tube simulation with those calcu- 
lated for the three substitute gas mixtures for an a + p of 4°. 
For comparison we have shown the pressure distribution that would 
be obtained by expanding air at a constant y of 1.4, assuming 
the exit Mach number matches that of the actual engine. 

We have integrated the pressure distributions shown in 
Figs. 13 through 15 to obtain relative values for the forces in 
the normal and axial direction. We normalize the two dimensional 
forces by P 3 A, where A is the (unit) width times L, the 
afterbody length. We have not calculated any pitching moment 
changes, because the shifts in center of pressure are extremely 
small. Figure 16 is a bar graph comparing all of the forces cal- 
culated for the case of a + p = 4° , and Fig. 17 is a comparison 
of those for the a + p = 12° case. The differences in normal 
force between frozen and equilibrium flow appear to be more sig- 
nificant than a quick look at the pressure distributions would 
indicate. This apparent disparity exists because the net normal 
forces are only 20 to 30 percent of the individual (cowl or 
afterbody) normal forces, and the percentage change in pressure 
level introduced by chemical reaction is much greater on the 
afterbody than it is on the cowl lip. The detonation tube simu- 
lation falls between the two cases. The axial forces are all 
within a few percent of one another with the exception of the 
constant y * 1.4 air case. Of the three substitute gas mix- 
tures chosen for this comparison, 50 percent Freon 13B1 and 50 
percent argon can be seen to give forces very close to the equi- 
librium full scale engine, and the detonation tube simulation 
lies between the frozen and equilibrium values. Once the true 
chemical state of the real engine exhaust flow field is known, 
there should be no trouble selecting a substitute gas to match 
the pressure distribution to any required system fidelity. 
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Fig. 15 Simulated Afterbody Pressure Distributions, 
M oo = 8, a + 0=4° 
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Relative Normal Forces 
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Fig. 16 Normal and Axial Forces for 2D Expansions 
M«o = 8 . or + (3 = 4° 
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Fig. 17 Normal and Axial Forces for 2D Expansions 
M *= 8, a + /3 = 12° 



CONCLUSIONS 


There are several substitute gas blends that should give 
excellent representations of the pressure, force, and moment 
distributions for scramjet exhaust flows over vehicle afterbodies 
These binary gas mixtures can be used in wind tunnels and in a 
shock tunnel, and they show excellent promise for representing 
the viscous behavior of the SCRJ flow as well as its inviscid 
characteristics. Different blends can be designed to compensate 
for small differences in the expected chemical reaction rates of 
the flight engine or for changes in flight conditions. 

The detonation tube system provides the best available lab- 
oratory simulation of the chemical, dynamic, and thermodynamic 
behavior of the flight SCRJ exhaust. This enables the detonation 
tube to be used as a standard of comparison against which substi- 
tute gas representations of the flight pressure distribution can 
be evaluated. 

Analysis shows that relatively small changes in flow thermo- 
chemistry and resultant pressure distributions can produce sig- 
nificant changes in vehicle loading, especially in afterbody 
normal force. This means that the efforts to overcome the dif- 
ficulties of creating and employing accurate methods of matching 
pressure distributions are, in fact, necessary, and that even 
relatively minor three dimensional flow details will need careful 
study in later stages of development. 

Proof-of-concept tests are required to show that the analyti 
cally predicted matching of the SCRJ exhaust can be realized in 
practice. These tests will also indicate the degree to which the 
detonation tube system can be used to investigate the effects of 
three dimensional flow details such as the dividers between seg- 
ments of the combustor. 

There is a need to develop a better understanding of the 
similitude requirements for the viscous mixing layer between in- 
ternal and external flow. If complete fidelity to all possible 
similitude requirements is enforced, the external flow will have 
to be of very low temperature or very high molecular weight, e.g. 
chilled Freon gas. It is likely that these constraints could be 
relaxed significantly to allow a much simpler test system, but 
only after experiments on the hypersonic mixing layer have shown 



this to be permissible. We suggest that another set of experi- 
ments be planned and conducted while the detonation tube proof- 
of-concept is under way, so that ultimate wind tunnel testing 
will be as simple as possible without loss of validity. 
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APPENDIX 


CALCULATION METHODS FOR SUBSTITUTE GAS MIXTURES 


A. Thermodynamic Properties 

The usual thermodynamic properties of homogeneous mixtures 
of two or more gases are computed from the component gas proper- 


ties 

by use 

of the relationships 

J 

a) 

ii 

S 

2X.c 
1 P i 


ratio of specific heats 

2X.(c - R) 

i p 


b) 

II 

2 X.ft. 
r i 

- 

molecular weight 

c) 

H 

m 

2X.H. 

l l 

- 

enthalpy 

d) 

S 

m 

2X.S. 

l l 

- 

entropy 

e) 

ii 

S 

a. 

0h B 

RT 

m m 

- 

gas law 


where X^ is the mole (volume) fraction of each component of the 
mixture . 


B. Transport Properties 

The transport properties, viscosity, thermal conductivity, 
and mutual diffusivity of binary and ternary mixtures are deter- 
mined by mixture estimation formulas. 

For viscosity, the Wilke method (Ref. A-l) has been employed 
because it yields results rapidly and has been verified by many 
experimental measurements and found to give comparable accuracy 
to the rigorous kinetic- theory relations (Ref. A-2) . The Wilke 
relation for a binary mixture viscosity is: 


1 

1 + X 2 9 21 


H- 


M- 

m 


X 2 

1 + X 1 q, 12 
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where 


(X = viscosity, (gm/cm - sec) x 10 (micropoise) 

= mole fraction of gas component 1 and 2, respectively 



m 


= molecular weight of gas component. 


Average deviations of computed values from data for many binary 
mixtures have been reported as typically one percent for the Wilke 
equation; this is considered adequate for our purposes. 


The viscosity estimation equation for a ternary mixture by 
Wilke's method takes the general form: 


M- 
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■ 1 tt i 
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1+5 C P- .(X./X.) 
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where and have been previously defined, and 
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We have selected the Mason-Saxena (Ref. A-3) estimation 
method for determining thermal conductivity, 1^, of binary mix- 
tures. This method is a simplified version of the classic Was- 
siljewa formulation, has the similar convenient form of the Wilke 
equation for mixture viscosity, and yields, typically, errors of 
less than 2 percent from actual data. The Mason-Saxena equa- 
tion for binary mixtures is: 



with units of cal/cm - sec - °K where the subscripts 1, 2, m 
refer to the two gas components and the binary mixture properties, 
respectively. 

Another transport property needed for two similarity parame- 
ters is the diffusion coefficient for a binary gas system. This 
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coefficient is independent of component gas concentration and 
quantifies the mutual diffusion of the gas mixture components for 
an isothermal process. 

The binary diffusion coefficient, D-^, ^ as been computed 
with the theoretical equation derived from modern nonuniform gas 
theory, using a Lennard-Jones model to characterize the mole- 
cules, where: 


D 12 = 0.001858 T 


3/2 [ C^i + 


Pa 12 fl D , 


2 , 

cm /sec 


where P is pressure in atmospheres and T is temperature in 
°K. The diffusion collision integral, n D , is tabulated 
(Ref. A-2) as a function of the normalized temperature, T/(e/k), 
where the molecular attraction law constant, e/k, for the mix- 
ture is obtained by use of the combining rule: 

(e/k) 12 =J (■e/k) 1 (e/k) 2 j 2 


Similarly, the molecular collision diameter constant, a-^, f° r 
the gas mixture is determinable from the rule: 

a 12 = ^ (a l + ° 2 ) 

where the component a's as well as (e/k)'s were obtained 
principally from Svehla's (Ref. A-4) listing. 

As a result of the computational investigation we have per- 
formed, the effect of the monomolecular gas additive in binary 
mixtures with Freon compounds can be characterized for several 
similarity parameters as follows: 

• For sufficiently high initial temperature to avoid 
phase change during an expansion process, increasing 
the additive increases the value of y. 

• For a given initial temperature, pressure, and model 
size, increasing the additive fraction decreases the 
Reynolds number. As a result the model size for match- 
ing the full scale Reynolds number of the SCRJ increases 
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as additive is increased at the same pressure. For- 
tunately, these blends are all much higher in molecular 
weight than air, so they produce higher Re at small 
scale, and matching is not difficult. 

• For a given initial temperature, the Mach number is 
matched at increasingly greater flow velocities in the 
model as the additive is increased. However, the model 
exhaust velocity is much lower than in the full scale 
apparatus. 

• For a given initial temperature, the Prandtl number 
decreases with increase in additive fraction. The 
substitute gas mixtures have a similar but slightly 
greater magnitude of Prandtl number compared to the 
SCRJ combustion products. 

• For a given initial temperature, the Schmidt number in- 
creases with increase in additive fraction. Both the 
substitute gases and SCRJ combustion products have 
similar magnitude of Schmidt numbers with the latter 
slightly greater. 

• The transport properties of viscosity and thermal con- 
ductivity increase with increase in additive fraction. 
The viscosity of a substitute gas in its proper 
operating range is generally about one third that of 
the SCRJ exhaust, and the thermal conductivity is more 
than an order of magnitude smaller than that of the 
hot SCRJ gases. 

• The mutual diffusivity of binary mixtures is greater 
than the single Freon-type component but is independent 
of mixture composition. The hot SCRJ exhaust gases 
have about two and one half orders of magnitude greater 
diffusivity than the relatively cool binary mixtures 
considered as substitute gases. 

C . One Dimensional Gas Dynamics 

The one dimensional gas dynamic properties of a substitute 
gas mixture flow have been determined with a minicomputer pro- 
gram, assuming an isentropic process and a thermally perfect gas. 
To justify these assumptions for the calorically imperfect (i.e., 
y 4 constant) binary gas mixtures, an over-all flow process has 
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been calculated on the basis of a series of small sequential 
step changes in temperature. An average y for the incremental 
change is determined, and used in calculating the new flow 
properties at the end of the step change, using the equations: 

P 2 /P l = ( t 2 / t i ) 7/Y ' 1 
1 

P 2 /Pi = (T 2 /T 1 )' Y " 1 

where 1 and 2 denote the beginning and end conditions, re- 
spectively, of the chosen small step interval in temperature: 


and 


7 


+ Y 2 


It is further assumed that across the small step, T 0 ^ - T o 2 ’ 
and the new Mach number at the end of the step change is found by: 



M 


m: 


i 

T 


The calculation starts with inputs of initial values of M, y , T, 
and p based on matching of substitute gas mixtures to SCRJ ex- 
haust gas properties. Generally, the Mach number is taken as the 
value at the combustor exit station, assuming homogeneous flow 
profiles. Initial values of y, T, and P are taken to match 
SCRJ combustor exit and Re, although the value of T can 

be selected for convenient test facility operation; the initial 
value of P can be taken as one and corrected for exact Re 
matching at a later time. Similarly, the value of initial y m 
can be taken as a value greater than the SCRJ y ^ so as to pro- 
duce an over-all substitute gas expansion process that better 
approximates the SCRJ process. 

The y m variation with temperature for a particular binary 
mixture is precalculated using Eq. (a) of Section A of this ap- 
pendix, and the JANAF table data of Cp for different temperatures 



of each mixture component. These discrete calculated values of 
Y m are curve fitted by a third order polynomial in T, stored 
in the minicomputer, and retrieved in the gas dynamic program to 
yield values of y m for each new temperature in the stepwise 
calculations through the flow channel. 

Additional inputs needed to commence the computation are 
the individual values of t to be used and the total number of 
prescribed repetitive steps. For directing the calculation to 
increasing values of M beyond the initial value, t is taken 
as slightly less than unity (e.g., 0.97). For decreasing 

values of M toward the stagnation point (M = 0) , a value 
of t slightly greater than 1 (e.g., 1.01) is selected. 

As the calculation approaches M = 0, an ongoing test for the 
next step being projected beyond M = 0 (i.e., to a virtual 

negative M) stops the calculation at the current step. If a 
finer grid size is desired to obtain nearly stagnation state 
conditions, the computation can be restarted at any small sub- 
sonic Mach number with a value for x closer to 1 than pre- 
viously used. It has been found that a final x = 1.0001 
usually is sufficient to obtain stagnation state conditions to 
five figure accuracy in temperature. 

Table A-l presents representative results of the one dimen- 
sional gas dynamic calculations for three good substitute gas bi 
nary mixtures. These results have been used in the two dimen- 
sional method of characteristics computations described else- 
where in this report in the Two Dimensional Afterbody Pressure 
and Force Prediction section. 



TABLE A-l 


REPRESENTATIVE ONE DIMENSIONAL GAS DYNAMIC PROPERTIES FOR THREE SUBSTITUTE GAS MIXTURES 


Substitute 
Gas Mixture 

M 

Mach 

No. 

Molecular 

Wt. 

y 

T 

(°K) 

P. 

(N/m 2 ) 

p 

(atm) 


( 






5 



1.3088 

94 

.42 

1.1795 

539.123 

2.558 

X 

10 3 

2.5247 


1.7205 



1.1853 

488.066 

1.286 

X 

io 3 

1.2693 


2.0287 



1.1911 

446 .26 

7.126 

X 

10 A 

0.7033 


2.36 



1.1994 

400.00 

3.607 

X 

10 4 

0.356 

50% Freon 13Bl+ ( 

2.9691 



1.2181 

324.00 

1.065 

X 

10 4 

0.1051 

507„ Argon 











3.5295 



1.2407 

262.44 

3.486 

X 

io 3 

0.0344 


4.0662 



1.2667 

212.576 

1.274 

X 

io 3 

0.01257 


4.5885 



1.2970 

172.1865 

5.168 

X 

10 2 

0.0051 


5.1113 


r 

1.3304 

139.471 

2.330 

X 

io 2 

0.0023 

A 

1.4519 

72 

324 

1.2054 

533.756 

1.891 

X 

io 5 

1.866 


1.7051 



1.2092 

497.875 

1.233 

X 

io 3 

1.217 


2.0224 

* 


1.2153 

450.725 

6.829 

X 

10 4 

0.674 


2.3314 



1.2232 

404.0 

3.678 

X 

10 4 

0.363 

407. Freon 12 + 

2.6537 



1.2328 

360.0 

1.978 

X 

io 4 

0.1952 

607» Argon < 











3.2078 



1.2531 

291.6 

6.758 

X 

io 3 

0.0667 


3.7313 



1.2788 

236.2 

2.530 

X 

io 3 

0.02497 


4.2441 



1.3079 

191.318 

1.044 

X 

io 3 

0.0103 


4.7479 



1.3423 

154.97 

4.762 

X 

io 2 

0.0047 

V. 

5.0046 


r 

1.3600 

139.47 

3.344 

X 

10 2 

0.0033 


1.5466 

103 

.15 

1.1786 

401.98 

1.284 

X 

10 6 

12.67 


2.299 



1.1951 

329.45 

2.634 

X 

10 3 

2.60 


2.365 



1.197 

322.96 

2.270 

X 

io 3 

2.24 


2.655 



1.2056 

295.30 

1.196 

X 

io 3 

1.18 

587. Freon 13B1 + 

2.9025 



1.2142 

272.70 

6.961 

X 

10 4 

0.687 

427. Argon < 








A 



3.173 



1.2244 

249.35 ' 

3.941 

X 

io 4 

0.389 


3.4099 



1.2338 

230.27 

2.452 

X 

10 4 

0.242 


3.821 



1.2508 

200.33 

1.150 

X 

io 4 

0.1135 


3.9977 



1.2594 

188.72 

8.542 

X 

10 3 

0.0843 

V. 

4.06 

\ 

r 

1.260 

185.0 

7.751 

X 

io 3 

0.0765 
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